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Abstract 

We present the data on depolarization of polarized neutron beam and second harmonic of mag- 
netization (M2) for Ndi_ x Ba x Mn03 (x = 0.23, the Curie temperature Tq ~ 124 K; x = 0.25, Tc ~ 
129 K) manganites. The depolarization starts to develop below T* ~ 147 K > Tc for both samples, 
being larger in x = 0.25 compound. This evidences the arising of a ferromagnetic (F) cluster phase 
below T* and a growth of its relative volume fraction with increasing doping concentration that 
agrees with the previously published results of M2 study. A characteristic size of the F clusters 
and their concentration are estimated combining the neutron depolarization and M2 data for x = 
0.25 manganite. 
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The interest in the study of the doped perovskite manganites is due to their unusual mag- 
netic and electronic properties, some aspects of which are not well understood. An important 
problem is the origination of an inhomogeneous magnetic state above Curie temperature, 
T C) in these compounds [1]. 

Ndi-^Ba^MnOs is a series of doped manganites which exhibits a transition from param- 
agnetic (P) to ferromagnetic (F) state for 0.2 < x < 0.35 {Tc ~ 120 K) [2]. For x > xjm ~ 
0.3, these compounds show metallic behavior below Tc, whereas for x < 0.3 they remain 
the insulators in the magnetic ordered phase. The arising of the inhomogeneous magnetic 
state above T c in process of a development of the second order transition was observed in 
the insulating Ndi.^Ba^MnOs (x = 0.23, 0.25) single crystals. It is characterized by the 
appearance of a magnetic phase with the strong nonlinear properties in the weak magnetic 
fields below T*(^ Tc + 20 K) [3-7]. The unconventional behavior was attributed to the F 
metallic regions which originate in the P matrix. The F clusters as well as the unconven- 
tional behavior above Tc were found in the other doped manganites [8-13]. This suggests 
that this phenomenon is an intrinsic property of these compounds. At the same time, the 
properties of this clustered state are not well elucidated, specifically a characteristic size of 
the F clusters is under question. 

The data on the magnetic, structural and transport properties for x = 0.23 and 0.25 NdBa 
manganites were published earlier [3-7]. According to the second harmonic of magnetization 
(M 2 ) measurements, the signal arising below T* from the F metallic regions indicated an 
increasing in a volume fraction of the F phase at T — > Tc, this fraction being larger in the 
x = 0.25 compound. 

In this paper the data on a depolarization of a transmitted polarized neutron beam are 
presented for both manganites. Also the additional data on the field dependence of the M 2 
at T =142.6 K for x = 0.25 crystal are reported. The combined data on the depolarization 
and on M 2 (H) for this compound are used in the quantitative analysis. It shows that the 
single-domain F regions, which above T c are in a close to superparamagnetic (SPM) regime, 
can reveal the strong nonlinear behavior in the weak magnetic fields. The analysis allows 
us to estimate the parameters of these F clusters (a magnetization, a characteristic size, a 
concentration, a magnetic anisotropy and a relaxation rate of the magnetization). 

The same Ndi-^Ba^MnOs (x = 0.25) single crystal as in Refs. (6,7) was employed at 
the M 2 measurements. It was grown by the flux melt technique, using the BaO-B203-BaF2 
ternary system as a solvent [14]. Structural single phase of the crystal was confirmed by 
the X-ray and neutron diffractions. The cation composition of the crystal was determined 
by X-ray fluorescent analysis. The average oxidation state of manganese determined by 
photometry [15] was found to be close to the value expected for the stoichiometric oxygen 
content. The powder samples used in the depolarization measurements were prepared by 
the traditional manner [2]. Structural single phase of the samples was controlled by X-ray 
diffraction, and the oxygen content was determined by a thermogravimetric analysis. 

The polarization measurements were performed using a small-angle polarized neutron 
setup "Vector" with a wavelength A = 9.2 A, AA/A = 0.25 (WWR-M reactor Gatchina) 
[16]- 

The additional measurements of the second harmonic of the longitudinal component of 
magnetization M 2 were performed in the parallel steady and alternating magnetic fields 
H + hsmut(h < 37 Oe, uj/2tt 15.7 MHz) in the temperature range T* > T > Tc- The 
ReM 2 and ImM 2 parts of M 2 were simultaneously recorded as the functions of H. This field 
was scanned symmetrically relative to the point if =0 for detecting a field hysteresis of the 
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signal. The amplitude of iJ-scan was 300 Oe. A condition M 2 oc h 2 was satisfied in the 
measurements. An installation and a method of separation of the M 2 -phase components 
have been described previously [17]. 

Fig. 1 shows the temperature dependence of the polarization for Ndi-^Ba^MnOs x = 
0.23, 0.25 polycrystalline manganites. It is seen that in both samples the depolarization 
appears just below T* « 147 K > T c (T c « 124 K for x =0.23 and 129 K for x =0.25), 
indicating the origination of the F regions, and it increases at decreasing temperature. The 
depolarization is larger for x = 0.25 sample. These results confirm the M 2 (T, H) data [4,6]. 

The M 2 response of the F clusters coexists with the critical contribution of the paramag- 
netic matrix, which also increases with decreasing temperature. Therefore, it is convenient 
for the analysis to choose the M 2 (H,T) curves of the x = 0.25 compound at T = 142.6 K 
(well above Tc), which are displayed in Fig. 2. A characteristic kink in the ReM 2 (if ) depen- 
dence, which is due to a competition of the paramagnetic critical contribution and the signal 
of the F regions, is clearly seen at this temperature. The lmM 2 (H) is mainly determined 
by the F clusters [6]. Additionally, the critical paramagnetic contribution is proportional 
to H (see below) and the demagnetization corrections are small here. The subsystem of F 
clusters can be considered as an ensemble of single domain F particles in nearly superpara- 
magnetic regime since the field hysteresis of M 2 (H) dependences is small [4,6]. Nevertheless, 
the presence of this small ff-hysteresis can be seen in lmM 2 (H) component (Fig. 2b) as an 
incomplete inversion symmetry of the curve relative to the point H = 0. 

Let us go to the analysis. The neutron polarization P after passing of the sample with 
the ferromagnetic regions of a small concentration is given by [18] 

P = P exp[-A/3( ln B/u) 2 RC 1/3 L], (1) 

where Pq is the initial polarization directed along the beam, 7„ is the gyromagnetic ratio 
of the neutron, v « 4.3-10 4 cm/s is the neutron velocity, B = 4irgfiB < S > /Vq is the 
induction in the F clusters, < S > is the value of the spin at temperature T, Vq is the 
volume per magnetic atom (58.6 A 3 [7]), R is the mean radius of the F regions, C is the 
relative concentration of the superparamagnetic phase and L 0.2 cm is the thickness of 
the sample. We have three unknown parameters B, R, C in (1) and need the additional M 2 
results to find them from experimental data. 

Going over the M 2 analysis, let us consider a possible magnetic anisotropy of the F regions. 
Our compound has the Pbnm space group with a relationship between the lattice parameters 
a > b cj y/2 [7] that suggests a dominating uniaxial anisotropy directed along the a-axis. 
We assume an easy axis character of the anisotropy. This assumption is supported by our 
preliminarily magnetic resonance measurements in the x = 0.3 single crystal (with a close 
to the x = 0.25 manganite crystal structure) where the resonance signal from the F regions 
is observed above T c . Such a signal is not detected at x — 0.25 [7] because of a small 
concentration of the F clusters. As a result, the Hamiltonian of the single-domain F region 
is given by 

H = -N[0 + K(ee,) 2 }, (2) 

where N is the number of the spins in the region, \x = gfiB < S >,g ~ 2 is the g - factor, H is 
the steady magnetic field, K is the effective uniaxial anisotropy per a magnetic atom which 
can include a shape anisotropy, e is the direction of the uniaxial anisotropy and = \Lj\x. 
The M 2 (H) data, which obtained for the magnetic fields directed along the different crystal 
axes, do not reveal an orientation dependence. It suggests a twin structure with a nearly 
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equal population of the structural domains. We will show below that the domains with easy 
axis give the main contribution to M 2 (H) for H directed along the crystallographic axes. In 
this case, the longitudinal response of the second order for single-domain F cluster can be 
written as [4]: 

M(Hn)/h 2 - wm iu dT{H)/dH 

Here 2x2{H) = d 2 M/dH 2 ,xi(H) = dM/dH, M is the magnetization of the single-domain 
region and T(H) is its relaxation rate. The magnetization of the domain is given by 

M(mH/T, KN/T) = mL(H/C 1} a), (4) 

where m = fiN, a = KN/T, C\ = T/m and 

1 

L(H/C 1 ,a) = C 1 dhiZ/dH, Z = J coth(x///d) exp(ax 2 )dx. (5) 

o 

Here Z is the partition function. In the first approximation, we perform a fit of the data by 
neglecting the second net dynamic term in Eq. (3) and ignoring an if -dependence of T. The 
M 2 of the sample, which is the sum of the response from all the F clusters M 2S = CM 2 and 
a contribution from the paramagnetic matrix M 2 c, can be written as 

M 2 (H, oo)/h 2 = (C> 2 + iC»)x2(H)/m + (C 3 + X%)H, (6) 



with 



2C ' 2 ~ \^ c {2u/v) 2 + r c ' 2 - Y c ' 2: (7) 



QC 3 = (1 - C) 



d 2 M c 



dH 2 



H=0 T C 



C> = ^C 3 , at r c »2cu. (8) 



Here a factor (1/3) in Eq. (7) takes into account the twinning structure of the crystal, 
Mq is the magnetization of the paramagnetic matrix in a weak field when M 2 q oc H [4], 
Tq is the relaxation rate of the uniform magnetization of the paramagnetic matrix. The 
factor C' 3 corresponds to ReM 2C at u = because (2u/T c ) 2 ~ 7.8 • KT 4 at T w 143 K [6]. 
One can perform a fit of the ReM 2 and ImM 2 component of the signal using Eqs. (4)- (8) 
where Ci, C 2 , C 3 , C 3 , C 3 and a are the fitting parameters. It allows one to find the following 
characteristics of the F regions: C/u — C < S >,Nfu = N < S >,KN and T which are 
determined by C 2 ,Ci,a and C 2 /C 2 , respectively. The value of < S > can be found by 
exploiting the data on the neutron depolarization (see below). As a result, we obtain all 
the needed parameters in this approach. A combined fit of the ReM 2 , ImM 2 and neutron 
depolarization data gives d w 42, a w 3, ss 4.23 • 10 4 ,C7 w 6.75 • 10" 5 , < S >^ 0.595 
as well as T/u pa 5. We have ignored above an if -dependence of T. Now we consider an 
effect of r(if). At further evaluation a variation of a makes the fit worse and we exploit 
below the found above a = 3 value. Using the numerical data for T(H) which is governed 
by thermoactivation over the barrier related to the magnetic anisotropy [19], one can obtain 
an approximate expression at a = 3 

2T(H) = r ^ P [-3(l + yi + e ? [-3(l-^l 
1 + 0.25H 2 + 1.5H 2 (1 - H 2 ) 
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where H = H/2H A ,H A = Kf/j, is the anisotropy field, Hj^ 1, dT(H)/dH > and 
T(H)/T(0) pa 8. At a following fit we use expression (3) for ReM 2 and ImM 2 with T(H) (Eq. 
(9)) in Eq. (6) instead of the first term. The fitting parameters are Cfu = C < S >, C±, T 
and Cg. The M 2 data give Cfu = C < S > and jV/# = N < S >. As a result, one can find the 
< 5 > from Eq.(l) assuming a spherical form of the regions (R = ((?>/Ati)NVq) 1 ^ oc (Nf it / < 
S >)V3 oc 1/ < S > x / 3 , = (C m j < S >) x / 3 oc 1/ < 5 > 1 /3, B*RCVs oc< 5 > 4 / 3 ). 

Let us consider the experimental results. Fig. 2 presents the M 2 (H) dependences with 
the fitting curves (h ~ 36.6 Oe). The parameters of the F regions are found to be C\ ~ 47 
Oe, N w 3.26 • 10 4 , < S >^ 0.59, i? « 77.4 A, C w 8.3 • lO" 5 , F4 ~ 165 Oe, K « 1.31 • 10~ 2 
K and T{0)/u = T e~ 3 /u w 6.33, rjj« 624 MHz. We used P ~ 0.938 and P w 0.93 
at 142.6 K (Fig. 1). Note that ReM 2 (if) (Eq. (3)) is determined mainly by the static 
amplitude X2(H) because: (i) (2cj/T(0)) 2 pa 0.1 and r(if) increases with increasing H; 
(ii) as the calculations show, the contribution of the second dynamical term in Eq. (3) 
to ReM 2 contains a small factor [u/T(H)] 2 < 0.02 in comparison with the first one. At 
the same time, the imaginary parts of both terms have the comparable magnitudes and 
the similar iJ-dependences. We have neglected above a contribution of the twins with the 
anisotropy axes directed perpendicular to the field which is characterized by X2±(H). A 
ratio r(H) = 2x 2 j_(iJ)/x 2 (iJ) increases with increasing H, remaining, however, small in an 
essential region of changing X2(H) : r(H) < 0.1 from H = up to a maximum of |x 2 (if)| at 
H^ c pa 40 Oe and r(H) pa 1 only at 120 Oe where |x 2 ±(i?)| reaches a maximum. However, 
at this field \x2{H)\ becomes small X2(120Oe)/x 2 (iJ^ c ) pa 0.15 and ReM 2 (if) is determined 
here by the large linear critical contribution. As a result, the including of \x2±(H)\ does 
not affect the parameters of the F regions. A frequency dependence of X2±(w,H) was 
not analyzed. It is known a linear susceptibility Xi±( UJ ) H — whose ^-dependence is 
determined by a temperature spread of the resonance frequencies of the magnetic moments 
in the field of anisotropy [20] so that lmxi±(io) ~ (u/2K)xi± for to < 2K (pa 613 MHz). 
One can expect that in the magnetic field lmx2±(w,H) ~ (u/K)x2±(H) at the small uj 
and H < 2H A (pa 330 Oe). Since lm X2 (u,H) ~ 2(u/T(H)) X2 (H) and T(H) < 2K at 
H < 120 Oe we find, as for ReM 2 (uj, H), that the contribution of the "orthogonal" twins 
can only affect the tail of lmM 2 (u,H) at H > 120 Oe where r(H) > 1. These corrections 
are unessential because the parameters of the F domains are mainly determined by a region 
of H < 120 Oe and the linear critical paramagnetic contribution starts to dominate above 
120 Oe. Note, the taking into account the X2±(w,H) contribution to the response could 
probably improve the fit of the tails of the experimental curves but this question is out the 
scope of this paper. 

Let us discuss the results. The found attempt frequency r pa 1.25 • 10 4 MHz is within 
the typical interval of the values (10 3 — 10 5 MHz). The magnetization of the F regions 
characterized by < S > is small in comparison with a saturation magnetization (< S > /S ~ 
0.3, S = 1.875 is the average spin value for the given composition which corresponds to the 
saturation magnetization). It seems to be the natural result for the temperature slightly 
below T* pa 146 K. In the Lao.gSro.iMnOs manganite, the magnetization of the F clusters 
increased with decreasing temperature and was also small just below T* [10]. In regard to 
the magnetic anisotropy, it is known only an upper limit of the single ion magnetocrystalline 
anisotropy in the bulk (JaS 2 with J A < 0.5 K) from the ESR measurements above T c [7]. 
For the F regions, K can be written as K — J a < S > 2 , leading to J4 pa 0.04 K in agreement 
with the estimation. In addition, according to our preliminarily ESR data on the x = 0.3 
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NdBa single crystal, the anisotropy field of the F regions is Ha ~ 250 Oe which is not far 
from Ha ~ 165 Oe in our compound. The anisotropy fields of the F clusters were reported 
for La 9 Sr 01 MnO3 [10] and La 88 Ba .i2MnO 3 [21]. The dominating uniaxial components, 
which correspond to the hard and easy axes anisotropies in the first and second manganite, 
respectively, were found to be about 2.4 kOe. The small value of Ha oc J a < S > in our case 
(Ha/(2A kOe) ~ 0.07) can be explained by a smallness of < S > and J a- The presented 
data on the anisotropy fields [10, 21] were obtained at temperatures of 50 K below T* 
that suggested a larger value of < S > (2-3 times) than that at temperature of the several 
degrees below T*. The constant J a is related to the distortions of the oxygen octahedron 
and decreases with doping increasing. This constant can be essentially smaller (7-5 times) 
near the border of the insulator to metal transition in our system than that in the two 
above mentioned manganites where doping is slightly above a border of an antiferromagnet 
to ferromagnet transition. 

We obtain the rather large value for the size of the F regions R/a pa 19.9, where 
a = (Vo) 1//3 ~ 3.884 A is the magnetic lattice constant. It is essential, therefore, to compare 
it with a critical radius of a single F domain Rc/a ~ 2 • (SJ/ud) 1 ^ 2 (J is the exchange inter- 
action and ujy) = An{gix B ) 2 /V is the dipolar energy) for a weakly anisotropic ferromagnet 
when 2J A /oo D <C 1 (uo D pa 0.75 K and 2J A /^ D pa 0.035 in our system) [22]. Using the mean 
- field expression for Tc (~ T* pa 150 K) of the F regions, we find J ~ 14 K and Rc/a ~ 
15. This result suggests that the F regions are close to the uniformly magnetized domains. 

As it was discussed in Refs. (4,6), temperature evolution of the signal of the F regions 
corresponded mainly to an increase in their volume at decreasing temperature from T* down 
to Tc- The matter is that a field position of the extreme (H^ 1 30 Oe, Fig. 2b) in ImM2(T, H), 
which is due to the F regions, remains temperature independent with decreasing temperature 
whereas the extreme value of ImM 2 (T, H^) increases sharply down to T C - This behavior 
reflects an increase in the concentration of the F clusters without the noticeable changes in 
their parameters. The same peculiarity of the ImM2(T, H^ 1 ) dependence with H^ 1 pa 22 
Oe was also observed in the x = 0.23 NdBa crystal [4]. The closeness of H^ 1 and a strong 
resemblance between the ReM 2 (T, H) functions in both systems suggests that the F regions 
are characterized by the close values of the parameters. The main difference between the 
compounds is the concentration of the F clusters which is about 10 times larger at x = 0.25 
than that at x = 0.23 [6]. 

Magnetization of the F regions is characterized by the small scale of the field nonlinearity 
C\ pa 47 Oe (Eq. 4). As a result, the M(T,H) measurements at H — 1 kOe performed 
above T c in the x = 0.25 NdBa compound reveal only the critical temperature dependence 
of the paramagnetic matrix [7] since the magnetization of the F clusters, being near to a 
saturation for the superparamagnetic regime, gives a temperature independent correction 
which is proportional to their rather small concentration. An extreme sensitivity of the 
clustered state to magnetic field was found in Lai-^Ba^MnOs (x = 0.27 and 0.3) [13]. The 
anomalous contribution to xi was suppressed by H ~ 150 Oe in the paramagnetic region 
(for T = 275-350 K, T c pa 245 K, x = 0.27 and T = 310-345 K, T c pa 310 K, x = 0.3). 
This behavior can be explained by a small characteristic if-scale for magnetization of the F 
clusters which has the comparable value with one in our system. 

In conclusion, the analysis of the data on the nonlinear longitudinal response and the 
depolarization of a polarized neutron beam for the Ndo.75Bao.25Mn03 manganite allows us 
to estimate the characteristic size of the F clusters (R ~ 77.4 A), which form in critical 
paramagnetic regime below T* pa 147 K and lead to a formation of heterogeneous magnetic 
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state above T c . 
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Figure captions 

Fig. 1. Temperature dependence of the depolarization of the polarized neutron beam for 
Ndi_a;Ba x Mn03 (x = 0.23, 0.25) powdered samples. Neutron wavelength is A = 9.2 A, 
AA/A = 0.25 and the thickness of the sample is L 0.2 cm. 

Fig. 2. Second harmonic of magnetization, M2(H), as a function of H for the Ndi-^Ba^MnOs 
single crystal. Panels show iJ-dependence of the phase components ((a) - KeM2(H), (b) - 
lmM2(H)) at T = 142.6 K and their fits. The fits are described in the text. 
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